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As a generic model describing quasi-one-dimensional Mott and Peierls insulators, we investigate 
the Holstein-Hubbard model for half-filled bands using numerical techniques. Combining Lanczos 
diagonalization with Chebyshev moment expansion we calculate exactly the photoemission and 
inverse photoemission spectra and use these to establish the phase diagram of the model. While 
polaronic features emerge only at strong electron-phonon couplings, pronounced phonon signatures, 
such as multi-quanta band states, can be found in the Mott insulating regime as well. In order to 
corroborate the Mott to Peierls transition scenario, we determine the spin and charge excitation 
gaps by a finite-size scaling analysis based on density-matrix renormalization group calculations. 
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The one-dimensional (ID) Holstcin-Hubbard model 
(HHM) has been used extensively to describe for novel 
low-dimensional materials, e.g., conjugated polymers, 
organic charge transfer salts or halogen-bridged tran- 
sition metal complexes^, and the associated mctal- 
insulatoii^i^ and insulator-insulator transitions^i^^ The 
HHM accounts for a tight-binding electron band, intra- 
site Coulomb repulsion between electrons of opposite 
spin, and a local coupling of the charge carriers to the 
phonon system: 
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Here, c]^ {Ci„) denote fermionic creation (annihilation) 
operators of spin-cr electrons (tr =t, |) on a ID lattice 
with N sites, = c\^c^^, and b\ (6^) are the corre- 
sponding bosonic operators for a dispersionsless optical 
phonon with frequency loq. 

The physics of the model is governed by the com- 
petition between electron itinerancy {oz W = At) on 
the one hand and electron-electron (cx u = U/At) and 
electron-phonon (oc A = Sp/'it; Sp is the polaron shift) 
interactions on the other hand, which both tend to im- 
mobilize the charge carriers. At least for the half-filled 
band case (^^ ^ riia = N^i = N), Mott insulator (MI) or 
Peierls insulator (PI) states are expected to be favored 
over the metallic state at temperature T = (see Fig. 1). 
The correlated MI shows pronounced spin-density-wave 
(SDW) fluctuations but has continuous symmetry. It 
therefore exhibits no long-range order in ID. In contrast 
the PI is characterized by dominant charge-density-wave 
(CDW) correlations and true long-range order because a 
discrete symmetry is broken. While the gaps to both spin 
(As) and charge (Ac) excitations are finite in the PI, the 
spin gap vanishes in the ID MI, which is related to spin 
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FIG. 1: Schematic phase diagram of the Ifolstein-Hubbard 
model at half-filling. 

charge separation. In a strict sense these results hold 
in the adiabatic limit {ujq = 0) for "J7-only" (Hubbard 
model) and "A-only" (Peierls model) parameters. At fi- 
nite phonon frequency and U = (Holstein model) a crit- 
ical electron-phonon (EP) coupling is required to set up 
the CDW phase characterized by alternating doubly oc- 
cupied and empty sitesi The concomitant gap formation 
and metal-insulator transition has recently been studied 
in the limit of infinite dimensions^ Depending on the adi- 
abaticity ratio a = ujo/t the PI represents a traditional 
band insulator (a <IC 1) or a bipolaronic insulator (a 3> 1, 
= Ep/ujo ^ l)i^ Although for the more general HHM 
the situation is much less clear, we expect that the fea- 
tures of the insulating phase will depend markedly on 
the ratio of Coulomb and EP interactions u/X, allowing 
for quantum phase transitions between insulating phases. 
Indeed, based on recent numerical results for the stag- 
gered spin- and charge-structure factors, it has been ar- 
gued that the HHM shows a crossover between Mott and 
Peierls insulating phases at m/A ~ liifiMore precisely, for 
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finite periodic cfiains, the MI-PI quantum phase transi- 
tion could be identified by a ground-state level crossing 
associated with a change in the parity eigenvalue 
Note that this scenario differs from the (weakly interact- 
ing) HHM with frozen phononsiS, where there is strong 
evidence in favor of two quantum critical points, as in 
the ionic Hubbard modeliii*i^ 

The aim of this work is to establish the physical picture 
developed to date for the interplay of spin, charge and 
lattice degrees of freedom in the ID Holstein Hubbard 
model. In particular we attempt to verify the proposed 
phase diagram by examining the single-particle excita- 
tions. For these purposes we employ Lanczos exact di- 
agonalization (ED)^^, kernel polynomiaU^, and density- 
matrix renormalization group (DMRG)^^'^^ methods to 
determine the (inverse) photoemission spectra as well as 
the spin and charge excitation gaps. These quasi-exact 
numerical techniques allow us to obtain reliable results 
for all interaction strengths with the full quantum dy- 
namics of phonons taken into account. Exact diagonal- 
izations are seriously limited in achievable system sizes 
N, but have the advantage that spectral quantities are 
easily accessible. Complementary, the DMRG algorithm 
yields specific eigenstates of large systems by implement- 
ing a renormalization scheme and an optimal truncation 
of the Hilbert space. Thus it permits for a finite-size anal- 
ysis of the ground-state energies in different particle and 
spin sectors, which is required to determine the behav- 
ior of the various excitation gaps in the thermodynamic 
limit. 

We begin by studying the spectral density of single- 
particle excitations associated with the injection of a 
spin-cr electron with wave number K, A'^^{uj) (inverse 
photoemission (IPE)), and the corresponding quantity 
for the emission of an electron, A~^^{uj) (photoemission 
(PE)), where 
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with c^^ = c\^^ and c^^ = Cj^^. {ipo^"^) is the ground 
state of the system with N^i electrons and liprn'''^^^^) are 
eigenstates of the {Nei ± l)-particle system. Eq^"'^ and 

Em''^^'^^ are the corresponding energies. Adding the 
spectral densities of (photo-) emission and absorption 
we obtain the spectral function Aj^£,(a;) = A^^{ijj) + 
A~^^{lij), which obeys various sum rules and allows for a 
connection to angle-resolved photoemission spectroscopy 
(ARPES). The simplest sum rule, /^A;f<,(w)du; = 1, 
reflects the normalization of AKai'^) but is not useful for 
ARPES since it involves both occupied and unoccupied 
states. 'Y^^ ^_^nF{uj)AK(y{^)d^ = i^-iK) (where np{uj) 
is the Fermi function) is more important, since it relates 
the ARPES intensity to the number of electrons in a mo- 
mentum state K: n{K) = ^^icl^^Cj^^). The ED results 
presented for A^^(uj) in the following were obtained for 



an eight-site system with periodic boundary conditions^ 
Let us first consider the MI regime. Figure 2 dis- 
plays the IPE and PE spectra for the HHM at the 
allowed wave numbers of our finite system: K = 
0, ±7r/4, ±7r/2, ±37r/4, and tt. To reliably monitor a 
possible band splitting induced by the Hubbard and EP 
couplings at half-filling it is necessary to guarantee that 
the Fermi momenta Kp ~ ±7r/2 are occupied, which is 
the case for N — Al [l integer, periodic boundary con- 
ditions). The most prominent feature we observe in the 
MI regime is the opening of a gap at K = ±7r/2, in- 
dicating massive charge excitations. A comparison with 
the results obtained for the pure Hubbard model classi- 
fies this gap as the Mott-Hubbard correlation gap. Its 
value A/i ~ 3.25 almost coincides with the optical gap 
Aopt we determined by evaluating the regular part of 
the optical conductivity for the same parameters. The 
dispersion of the lower (upper) Hubbard band can be de- 
rived tracing the uppermost (lowest) excitations in each 
K sector. Due to the finiteness of our system and the 
rather moderate value u = 1.5, PE (IPE) excitations 
with K = ±37r/4 and tt (K ^ ±tt/4: and 0) have stiU fi- 
nite spectral weight. This can be seen from the integrated 
spectral densities S^^{uj) = J'^^dto' A^^{uj'), which, in 
addition to the sum rule >5'^^(— oo) + S^^{oo) — 1, sat- 
isfy the relations S'|^(±oo)-H5'^_^^^(±cx)) = 1 {K > 0). 
Since the spectral weight of the PE excitations with 
K > tt/2 is expected to vanish as N goes to infinity for 
M ^ 1, the lower Hubbard band will be completely filled 
CL\K\<Kf.,a 1^°^ d.ujA-^^{u!) ~ Nei), and consequently 
the system behaves as an insulator at T = OiS£ As a 
result of the coupling to the phonon system the elec- 
tronic levels in each K sector split, creating phonon side 
bands. The distinct peaks are separated by multiples of 
the bare phonon frequency and can be assigned to re- 
laxation processes of the Q = phonon modesiSi The 
number of phonons involved is controlled by g^. S^^{uj) 
shows clearly that the total spectral weight of the result- 
ing excitation bands equals the weight of the respective 
electronic excitations in the pure Hubbard model. In- 
terestingly, mediated by Q 7^ phonons, there appear 
"shadows" of the band belonging to a dominant elec- 
tronic excitation in a certain K sector in other K sectors, 
giving rise to a weak "breather-like" excitationSSi, which 
is almost dispersionsless in the Brillouin zone. 

If we decrease the Hubbard interaction at fixed EP 
coupling strength the Mott-Hubbard gap weakens and fi- 
nally closes at about {u/X)c — 1, which marks the MI-PI 
crossover. This is the situation shown in Fig. 3. Ap- 
proaching the critical point from above and below, the 
ground state and the first excited state become degener- 
ate. These states have different eigenvalues P of the site- 
inversion operator Pcj^P'^ = cl^_i„ {i = 0, . . . , N — 1) 
and we have verified that the ground-state site parity is 
P = -M in the MI and P = -1 in the PI. Obviously the 
critical point is characterized by gapless charge excita- 
tions at the Fermi momenta but should not be considered 
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FIG. 2: Wave- number-resolved spectral densities for pho- 
toemission {A]^^{uj); red lines) and inverse photoemission 
(A^^(tj); black lines) in the Mott insulating state {u/\ S> 1). 
The corresponding integrated densities S'l^^{u)) are given by 
dashed lines. Data for the pure Hubbard model (blue and 
green lines) were shifted by —{£pN^i/N) and included for com- 
parison. 




FIG. 3: PE (red lines) and IPE (red lines) spectra near the 
Mott insulator Peierls insulator transition point (u ~ A). 

as metallic because the Drudc weight is ill-defincdiiiiS^ 

If the Hubbard interaction is further reduced, i.e. the 
EP coupling overcomes the on-site Coulomb repulsion, a 
CDW accompanied by a dimerization of the lattice de- 
velops. As a result the electronic band structure be- 
comes gapped again (see Fig. 4 (upper panel)). The 
form of the spectra, however, is quite different from MI 
case. While in the MI regime the lowest peak in each 
K sector is clearly the dominant one, in the BI phase 
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FIG. 4: PE and IPE spectra in the Peierls phase (u < A). 
The upper (lower) panels show typical results obtained for 
the case of a band insulator (BI) at a <C 1 and bipolaronic 
insulator (BPI) at a ^ 1, respectively. 



rather broad (I)PE signatures appear. Within these ex- 
citation bands the spectral weight is almost uniformly 
distributed, which is a clear signature of multi-phonon 
absorption and emission processes that accompany every 
single-particle excitations in the PI. The lineshape then 
reflects the (Poisson-like) distribution of the phonons in 
the ground state. Again low-intensity "shadow bands" 
become visible. Remarkably, now the cumulative spec- 
tral weight J2\k\<tt/2 ^Kai~°^) gives nearly the total 
number Nei of electrons, i.e., the finite-size effect men- 
tioned above for the MI are much less pronounced in 
the BI state. The situation changes radically if the in- 
sulating behavior is associated with localized bipolarons 
forming a CDW state (see Fig. 4, lower panel). Due to 
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using DMRGii^ Here E^\s^) is the ground-state en- 
ergy of the HHM with M partcles in the sector with total 
spin- 2; component 5^. 

Obviously the finite-size scaling presented in Fig. 5 
for Af,/s substantiates our introductory discussion of the 
phase diagram (cf. Fig. 1). Ac and are finite in the 
PI and will converge further as TV ^ 00. Both gaps seem 
to vanish at the quantum phase transition point of the 
HHM with finite-frequency phonons, but in the critical 
region the finite-size scaling is extremely delicate. In the 
MI we found a finite charge excitation gap, which in the 
limit u/A ^ 1 scales to the optical gap of the Hubbard 
model, whereas the extrapolated spin gap remains zero^ 



FIG. 5: DMRG finite-size scaling of spin- and charge ex- 
citation gaps in the HHM with dynamical phonons (A = 
0.35, Q = 0.1). Note the different scale of Ac in the MI 
phase. Open and filled symbols denote DMRG results for 
periodic (PBC) and open (OPC) boundary conditions, re- 
spectively. The accessible system sizes are smaller at larger 
A/u, where an increasing number of (phononic) pseudosites is 
required to reach convergence with respect to the phonons. 
Stars represent the ED results for the eight-site system. 
The arrow marks the value of the optical gap Aopt for the 
Bethe-ansatz solvable ID Hubbard model, which is given by 



Aopt/4i : 



1 -|- ln(2)/2M in the limit of large u> Ir 



strong polaronic effects an almost flat band dispersion 
results with exponentially small (electronic) quasiparti- 
cle weight. Now the dominant peaks in the incoherent 
part of the (I)PE spectra are related to multiples of the 
(large) bare phonon frequency broadened by electronic 
excitations. 

Since many-body gaps to excited states form the basis 
for making contact with experimentally measurable exci- 
tation gaps and can also be used to characterize different 
phases of the HHM, we finally determine the charge and 
spin gaps. 
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In summary, we have presented a comprehensive pic- 
ture of the physical properties of the ID half-filled finite- 
phonon frequency Holstein-Hubbard model. With re- 
spect to the metal the electron-electron coupling favors 
the Mott insulating state whereas the electron-phonon 
interaction is responsible for the Peierls insulator to oc- 
cur. The PI typifies a band insulator in the adiabatic 
weak-to-intermediate coupling range or a bipolaronic in- 
sulator for non-to-antiadiabatic strong-coupling. Our re- 
sults for the single-particle spectra and spin/charge ex- 
citation gaps give clear indication of a Mott- to Peierls- 
insulator quantum phase transition at u/\ ~ 1. Quan- 
tum phonon dynamics yields pronounced effects in the 
(I)PE spectra, which might be of great importance for in- 
terpreting photoemission experiments of low-dimensional 
strongly correlated electron-phonon systems such as MX- 
chain compounds^ 
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